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Abstract: The equilibrium properties of a tangent-hard- sphere polymer chain in a hard 
sphere monomer solvent are studied using a Born-Green- Yvon integral equation in 
conjunction with a two-site solvation potential. The solvation potential is constructed using 
low density results for a hard-sphere trimer in a hard sphere solvent. The BGY equation 
has been solved for polymers of lengths up to 100 for a range of solvent densities and 
solvent diameters. The theory accurately describes the compression of the average polymer 
dimensions with increasing solvent density. The theory also accurately describes the 
variation in this compression as a function of the solvent diameter, predicting a maximum 
effect when the solvent diameter is two to three times larger than the polymer site diameter. 
Scaling exponents relating the polymer dimensions to chain length and solvent density are 
also obtained. 



I. INTRODUCTION 



The statistical properties of a polymer in solution are strongly dependent on the nature 
of the solvent [1,2]. In a so called "good" solvent the dimensions of the polymer are 
expanded relative to those of a non-interacting or ideal polymer while in a "poor" solvent 
the polymer dimensions are collapsed relative to the ideal or theta state. A flexible hard- 
sphere chain with site diameter a in a fluid of hard sphere monomers with diameter D 
provides a model for a polymer in a solvent. Even in this simple athermal system 
interesting solvent effects are observed [3-6]. For D/a<l the hard sphere fluid represents a 
good solvent for the hard-sphere chain at all densities. However, for larger D/a Suen et al. 
[6] have recently reported evidence that a high density hard sphere fluid provides a poor 
solvent for a long hard- sphere chain. In any case, with increasing solvent density the 
average dimensions of the hard-sphere polymer are compressed. Furthermore, the 
magnitude of this compression effect is strongly dependent on the ratio D/a. 

One approach for studying a n-mev chain in solution is to introduce a rc-site solvation 
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potential which accounts for the interactions of the solvent with the polymer [5,7,8]. 
Approximating the n-site potential as a sum of two-site potentials allows one to map the 
polymer- solvent problem onto that of an isolated polymer whose segments interact via an 
density-dependent effective potential. In this work we introduce a two- site solvation 
potential and study the effects of solvent size and density on a hard- sphere-chain polymer 
using a Born-Green- Yvon (BGY) integral equation [9]. 



H. THEORY 

A. Intramolecular Distribution Functions 

In this work we study the equilibrium structure of a single polymer molecule in a 
solvent. The polymer is modeled as a flexible chain of n identical hard- sphere monomers 
with diameter a which are connected by universal joints of bond length a. The solvent 
consists of N hard spheres of diameter D and the solvent density is p=N/V where V is the 
total system volume. The configurational properties of the polymer can be expressed in 
terms of a set of intramolecular site-site distribution functions w y <r; p,D) which are related 
to the probability that two sites on the polymer, i andj, are separated by a distance 
r= \ x-v j \. The exact expressions for the w ; y functions are quite formidable since they 

involve integrals over the complete set of the N solvent molecule coordinates. We note that 
this full problem is amenable to a BGY type integral equation treatment [9,10] although 
here we take a much simpler tack. 

A formal simplification in the definition of the Wy functions can be achieved by writing 
the integrals over the solvent coordinates in terms of an n-site solvation potential [7]. Now, 
by making the assumption that the rc-site solvation potential can be written as a sum of two- 
site solvation potentials (Chandler and Pratts' superposition approximation) we can map the 
full problem described above on to that of an isolated chain with an modified site-site 
potential which implicitly accounts for the effects of the solvent. The resulting density 
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dependent effective potential between nonbonded polymer sites i and j (\i-j\>\) is given by 
uf (r; p,D) = u)f(r) + u s f(r; p,D) (1) 

where u^f(r) is the original hard sphere potential 

y 



= j _ (2) 
y [ r>a 

and u s ^(r; p,D) is the site-site solvation potential. 

The site-site distribution function for an n-mer chain whose segments interact via such 

an effective potential is given by 

Wy(r; p ' D) = vT\ f f n n ex p[-K?(r; p'°)] n wi n ^ (3) 

where (3=l/k B T, s(r) = 6(r-a)/4jta 2 is the intramolecular distribution function between 
bonded sites, and Z n (p) is the single chain partition function, which ensures the 
normalization condition 

J dr 4jtr 2 w^(r; p,D) = 1 . (4) 

Due to the constraints of chain connectivity the w -(r; p,D) functions are identically zero for 

y 

r>\i-j\o. A number of equilibrium configurational properties of the polymer chain may be 
expressed using averages over these distribution functions. In particular, the mean-square 
distance between sites / and j is 

(r.. 2 ) = j dr r 2 4jt r 2 w.<r; p,D) (5) 

and the mean- square radius of gyration is 



oo r<a 



n 
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B. BGY Integral Equation 

As we have shown previously [9], an integral equation for the intramolecular site-site 
distribution functions of an isolated interaction- site polymer can be obtained following the 
method of Born, Green, and Yvon [11,12]. The BGY equation for the end-to-end 
distribution function of an n-mer chain is given by 

%w ln = VjH3u^ f (r; p,D)] w ln + Jd7 2 VjlnsJ w^ 

n-l 

+ 2 fa Vj [-|3u^(r; p,D)] w (7) 
v=3 J 

where w lvn is a three-site distribution function which satisfies the reduction condition 

w, = fdr w, . (8) 

In J v Ivn v ' 

In order to solve the above BGY equation we must introduce two approximations. First, 
the required three-site functions are expressed in terms of two-site functions using a 
superposition-like approximation and, second, the distribution function between any two 
sites / and j is approximated by the end-to-end distribution function of a chain of length 
i+j-l. With these approximations, Eq. (7) can be solved recursively by direct numerical 
integration for a chain of arbitrary length. For details of the numerical solution the reader is 
referred to Ref. 9. 

This BGY theory has been used to study both isolated hard-sphere and square-well 
chain and ring polymers [9,13]. The success of the theory in these cases provides 
validation of the above approximations. In this work we apply the BGY theory to study a 
polymer in solution via the solvation potential approach described in the previous section. 
The construction of the required solvation potential is described forthwith. 
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C. The Solvation Potential 

Here we construct a polymer solvation potential which is motivated by low density 
results for a hard-sphere 3-mer in a hard sphere solvent. This task is accomplished by 
equating a low density expression for the hard-sphere 3-mer end-to-end distribution 
function to the analogous distribution function for an isolated 3-mer whose sites interact via 
an effective potential. 

To first order in solvent density, the end-to-end distribution function for a hard-sphere 
3-mer is given by 

exp|_Bu hs (r)l 

w 13 (r;p,D) = 13 [wg)(r) + p wfj^r; D)] a < r < 2a (9) 

z 3 (p) 

where the zeroth order coefficient is 

(0), 



W 13 



(r) = f dr 9 S.9S90 = — - — a<r<2a (10) 

J 2 12 23 8jta 2 r 



and we write the first order coefficient as 

w/Po-; D) = A 3 wg>(r) [5 3 (r; D) + 6 4 (a,a,r; D)]. (11) 
The /-bond diagrams comprising this first order coefficient are defined as 



and 



6 3 (r 13 ; D) = j d7 4 f u f M = ^ (16X 3 - 12rX 2 + r 3 ) a < r < 2X (12) 

6 4 (r 12 ,r 23 ,r 13 ; D) = J d7 4 / 14 f 24 / 34 (13) 



where X=(a+D)/2 and the solvent-polymer site interaction Mayer /-function is given by 
r -1 r<X 

/(r)= . v (14) 

[ r>X . 

An analytic expression for 5 4 is provided in Ref. 14. While rigorously the scaling factor A 
appearing in Eq. (11) should be unity, here we define it as the ratio of the solvent diameter 
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D to the interaction length X (i.e., A=D/X). In Table 1 we compare results for the mean- 
square end-to-end distance of a hard-sphere 3-mer in a D=a hard sphere solvent computed 
using the Eq. (9) first order expansion with those from Monte Carlo (MC) simulation [5]. 
The expansion correctly yields the solvent induced compression with increasing solvent 
volume fraction r| = jtpD 3 /6, although the magnitude of this effect is underestimated. 

The end-to-end distribution function of an isolated 3-mer whose sites interact via an 
effective potential is given exactly [15] by 

exp[-(3u^f(r; p,D)] fm 
w n (r;p,D) = 13 wfffr) . a < r < 2a (15) 

Z 3 (p) 13 

Using the definition of u eff (r) given in Eq. (1) and equating the above two expressions for 
w 13 (r) [Eqs. (9) and (15)] yields the following solvation potential 

(3u sol (r; p,D) = -In [ 1 + p wj^r; D)/w g>(r)] . a < r < 2a (16) 

In Fig. 1 we plot this site- site solvation potential both as a function of solvent volume 
fraction r) and solvent diameter D. In all cases the potential is short ranged and purely 
attractive. For a given solvent diameter the strength of the potential increases with 
increasing solvent density. The range of the potential decreases with decreasing D and the 
potential vanishes in both the limits of D— >0 and D— >oo. In the following we apply this 
solvation potential to the problem of a hard- sphere polymer in a hard sphere solvent. This 
potential mimics a very local, solvent induced site-site interaction mediated by a single 
solvent molecule. The potential to some extent incorporates intramolecular screening 
effects and contains information about the solvent geometry. The potential does not 
include any information about the local solvent structure. 
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III. RESULTS 



Using the single chain BGY equation [Eq. (7)] and the above polymer solvation 
potential [Eq. (16)] we have computed the intramolecular distribution functions of hard 
sphere chains in a hard sphere solvent for chains of lengths 4 < n < 100, solvent diameters 
0.1 < D/a < 100, and solvent volume fractions 0.0 < r| < 0.5r|(=jtpD 3 /6). Rather than 
show the actual w^(r; p,D) functions, here we present averages over these distribution 

functions in the form of the polymer mean-square end-to-end distance and radius of 
gyration defined in Eqs. (5) and (6), respectively. 

In Fig. 2 we show the variation in average polymer size with solvent density. BGY 
results for the mean-square end-to-end distance are shown for chains of lengths n=!0, 20, 
and 30 and solvent diameter D=a. Also shown in this figure are MC data for the polymer 
solvent system [4,5] and, in the case of n=20, MC data for an analogous polymer melt 
system [16]. The BGY theory is seen to correctly predict that the polymers are moderately 
compressed with increasing solvent density. For isolated chains (i.e., r|=0.0) the BGY 
theory slightly overestimates (r ln 2 ) and this discrepancy is seen to persist with increasing 

solvent density. Overall the BGY theory gives quite good agreement with the MC data 
even up to very high volume fractions. Note from Fig. 2 that both the average size, and its 
variation with density, of a hard-sphere 20-mer in a melt of hard-sphere 20-mers (filled 
symbols) is quite similar to the corresponding results for a hard-sphere 20-mer in a solvent 
of hard- sphere 1-mers. 

For a long chain in a solvent one expects a scaling relationship of the form R 2 ~n 2v 
(where R is either the mean- square end-to-end distance or radius of gyration), where an 
exponent of 2v = 1 corresponds to a theta solvent while 2v > 1 and 3/ < 1 correspond to 
good and poor solvents, respectively [1,2]. We have computed this exponent for hard- 
sphere chains in a hard sphere solvent with D=a. Our BGY results for (R g 2 ), for 
20 < n < 100, give exponents of 2v = 1.27, 1.21, and 1.16 for solvent volume fractions of 
r\ = 0.0, 0.3, and 0.5, respectively. The fact that the rpO.O result exceeds the Flory value of 
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6/5 indicates that the chain lengths considered here are not large enough to be in the 
asymptotic scaling regime. However, we can still conclude that the D=a solvent is a 
"good" solvent for all densities and that increasing solvent density corresponds to making 
the solvent conditions more "theta-like". This latter finding suggests that an athermal 
1-mer solvent produces similar effects on the dimensions of a single n-mer chain as an 
athermal n-mer solvent (for which 2v-»l with increasing system density [17]). 

For long chains in a dense solvent one also expects the scaling relationship R ~r\~^ with 
7=0.25 for a good solvent [18]. Although this value of y presumes that one is in the ideal 
scaling regime (i.e., 2v=l), we have nevertheless computed it for our hard-sphere 
chain/solvent system. Our BGY results for a hard-sphere n-mer in a hard sphere solvent, 
with D=a and 0.3<r|<0.5, give exponents of 7=0.3 and 0.4 for n=50 and 100, 
respectively. The MC melt data, shown in Fig. 2, give a scaling exponent of y=0.24 
(r)>0.2) while the BGY results for the polymer solvent system with n=20 give y=0.26. 

In Fig. 3 we show both BGY and MC results for the mean-square end-to-end distance 
as a function of solvent diameter D (0.1 < D/a < 100), for a hard-sphere 30-mer in a hard 
sphere solvent at solvent volume fraction r\ = 0.2. The BGY theory predicts that in the 
limits of a very large and very small solvent diameter D/a the average polymer size is 
essentially unperturbed from the isolated chain result. (We note that the average size of an 
isolated hard-sphere 30-mer is approximately equivalent to a sphere of diameter 7a). 
When the solvent diameter D is similar to the polymer site diameter a the polymer is 
compressed by the solvent with the maximum effect occurring for D/a ~ 2 to 3. In 
comparison with the Monte Carlo results [4], the BGY theory is found to be quite accurate. 



IV. CONCLUSION 

In this work we have studied the equilibrium properties of a hard- sphere polymer chain 
in a hard sphere monomer solvent. The polymer- solvent system is treated in the context of 
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a polymer site-site solvation potential. Using a BGY integral equation we have computed 
the intramolecular distribution functions for polymers over a range of solvent diameters and 
densities. The BGY theory accurately describes the compression of the average polymer 
dimensions with increasing solvent density. Additionally, the BGY theory accurately 
describes the variation in this compression as a function of solvent diameter, predicting a 
maximum effect when the solvent diameter is two to three times larger than the polymer 
site diameter 
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Table 1. Mean-square end-to-end distance for a tangent hard sphere 3-mer in a hard sphere 
solvent (D=a) computed via a first order density expansion and via Monte Carlo simulation 
[5]. 
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Figure Legends 



Fig. 1. Solvation potential as a function of site-site separation r over a range (a) of solvent 
volume fractions r\ with D=a and (b) of solvent diameters D with rp0.4. 

Fig. 2. Variation of the mean- square end-to-end distance with solvent density r| for a hard- 
sphere n-mer chain in a hard sphere solvent of diameter D=a. Solid lines are the results of 
the BGY theory and the open symbols are MC data from Gracye (() ) [5] and from 
Escobedo and de Pablo (O, □ ) [4]. The filled symbols are MC data for an n=20 polymer 
melt system taken from Yethiraj and Hall [16]. 

Fig. 3 Variation of the mean-square end-to-end distance with solvent diameter D for a 
hard-sphere 30-mer chain in a hard sphere solvent of volume fraction r\=0.2. Solid lines 
are the results of the BGY theory and the symbols are MC data from Escobedo and de 
Pablo [4]. The dashed line indicates the BGY result for an isolated hard- sphere 30-mer. 
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